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Abstract: The well-known M4 processes of Smith and Weissman are very flexible models for asymptotically 
dependent multivariate data. Extended M4 of Heffernan et al. allows to also account for asymptotic indepen- 
, dence. In this paper we introduce a more general multivariate model comprising asymptotic dependence and 

(N- independence, which has the extended M4 class as a particular case. We study properties of the proposed model. 

^ , In particular, we compute the multivariate extremal index, tail dependence and extremal coefficients. 
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^ I 1 Introduction 

I Smith and Weissman (|19j. 1996) presented the so called multivariate maxima of moving maxima (hence- 

■ forth M4) process for the modeling of cross-sectional and serial extreme dependencies. The class of 

M4 processes is very flexible for data which exhibit asymptotic dependence. However, there are also 
many examples of data which are asymptotically independent (see Ledford and Tawn [101 [TT] 1996/1997, 
among others), and hence cannot suitably be modeled by M4. Extended M4 processes (EM4) were latter 
^ ■ considered in Heffernan et al. ([5], 2007) in order to also account for asymptotic independence. 

\ The concept of tail independence between two random variables (r.v.'s), Xi and X2, with identical 

CO ■ marginal distribution function (d.f.) F, was introduced in Sibuya ([TT], 1960). More precisely, they are 

\l I said to be asymptotically independent if 

O ■ A = lim P{X2 > x\Xi > x) (1) 



is null and asymptotically dependent whenever A > 0, with xp ~ sup{a; : F{x) < 1} being the upper 
end-point. Coefficient A quantifies the amount of dependence of the bivariate upper tails and is usually 
denoted tail dependence coefficient (TDC). For instance, gaussian random pairs have A = 0, i.e., are 
■ asymptotically independent whereas i-distributed ones have A > and thus asymptotically dependent. 

?-H I In the asymptotic independent case, Ledford and Tawn ([TOl E] 1996/1997) proposed to model the 

null limit in ([T]) by introducing a new coefficient [t]) to rule the decay rate of the joint bivariate survival 
function evaluated at the same large x. More precisely, 

P{Xi >x,X2 > x) r-. >2;)^/'', asx^XF, (2) 



PiXi > x) 

where L is a slowly varying function, i.e. L(tx) / L{x) — > 1 as a; — > 00 for any fixed t > 0, and 77 G (0, 1] is 
a constant. The r.v.'s Xi and X2 are called positively associated when 1/2 < 77 < 1, nearly independent 
when 77 = 1/2 and negatively associated when < 77 < 1/2. Observe that they are asymptotically 
dependent if 7; = 1 and L{x) t^- 0, as x — > 00, and asymptotically independent otherwise. Explicit 
formulas for 77 can be seen in Heffernan ([7], 2000) for several known joint distributions. 

The lag-r (r G No) tail dependence coefficient for stationary d-dimensional sequences, {X„ = (X„.i, X„.(;)}n>i, 
with identical marginal distribution is naturally stated as (see Heffernan et al. 2007 and references 
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therein) 

A^-';)(X) = lim P{Xi+r r > x\Xi , > x). (3) 
Analogously, the lag-?- (r G Nq) Ledford and Tawn coefficient rjj^j}(X.) is defined by 

> x,Xi+r,f >x)^ l(^—-L__.^p(^Xi,, > x)'/''n'^^\ asx^ xp- (4) 

We have that A^y (X) = (X) is the TDC between the jth and the j'ih components, A^''^ (X) = A^'^'' (X) 

is the lag-r TDC within the jth sequence and A^y (X) is the lag-r cross-sectional TDC between the jth 
and the j'th sequences. A similar deduction concerns the Ledford and Tawn coefficients, respectively, 
,7,y(X),r;M(X) and ryg(X). 

A phenomenon also noticed in real data is that extreme events often tend to occur in clusters. The 
measure able to capture the clustered extremal dependence is the so called extremal index (Leadbetter 
et al. [S] f983). We shall define the multivariate extremal index (Nandagopalan jT2], f990) from which 
we can derive the univariate one of each marginal component. If {X„ = {X„.i, Xn_d)}n>i is an i.i.d. 
sequence such that, for some sequences of constants, {a„ = {a„i > 0,...,and > 0)}„>i and {b„ = 
(&„i, bnd)}n>i, the vector of componentwise maxima M„ = (-A/„ i, M„d) satisfies 

i^2(a„x + b„) = P(M„ < a„x + b„) H{xi, ...,Xd), (5) 

with H a non degenerate d.f., then iJ is a multivariate extreme value d.f. (MEV) and we say that 
F-g belongs to the domain of attraction of H, in short £ 'D{H). We recall that the convergence 
in ^ can also be stated using the copula concept, i.e., Cp^iFj^^^Xi), Fj^^{Xd)) = F^{xi, Xd), 
{xi,...,Xd) eR"^. It holds F^ e V{H) if and only if F^^ e T^{Hj), j = 1, d, and 

Cp;. (u^", ...,ity") — > CH{ui,...,Ud). 

A stationary sequence {X„}„>i, having common distribution i^x = F^, has extremal index 0{t) = 

e{Ti,...,Td) e [0,1] when, for each t = (Ti,...,rd) 6 there exists {ul^' = ("i^^l <J < d)}n>i, 
satisfying 



"(l-^x„,(<^^)) — > T„j = l,...,d, (6) 



P(M„ < u^;)) 7(t) and P(M„ < u^) -/M'^^^ 

Vectors ul^^ satisfying ^ are usually denoted normalized levels. Just as in one dimension, the extremal 
index is a key parameter relating the extreme value properties of a stationary sequence {X„}„>i to those 
of the i.i.d. associated sequence {X„}„>i. If ([S]) holds and {X„}„>i has multivariate extremal index 
9{t), then 

P(M„ <a„x + b„) — > G{xi,...,Xd), 
and the MEV d.f. in the limit satisfies: 

G(xi,...,xd) (7) 

and 



G,{x,)=H^^{x,), 
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with 

'^ji^j) = -~^ogHj{xj), j = and 9j = \im 0{ti, ...,Td). 

— >0 

The multivariate extremal index, although dependent of r, satisfies the following property: 

%Tl,...,CTd) =0(Tl,...,Trf), VC> 0. (8) 

In a univariate context, 0^ is called the extremal index of the sequence {Xnj}n>i- The value of l/9j can 
be interpreted as the limiting mean number of exceedances of a threshold per independent cluster as the 
threshold increases. A unit 9j means no serial clustering and is a form of asymptotic independence of 
extremes. 

Here we propose the pMAX model constructed from a d-dimensional stationary sequence {X„}„>i = 
X with multivariate extremal index and a sequence of independent and identically distributed (i.i.d) ran- 
dom vectors {Z„}„>i = Z, whose marginals are transformed through a positive exponent. The variation 
of the exponent values allows to obtain models ranging from serial asymptotic independence to asymp- 
totic dependence. The pMAX sequence may also exhibit clustering of large values. In computing the 
multivariate extremal index, tail dependence coefficients and the extremal coefficient we shall encounter 
a very rich and wide class that includes the EM4 processes (Section [5]) . We end with some brief notes on 
the estimation within the new model (Section [3]). 



2 The pMAX model 

Consider {X„ ~ (^ri,i, ^n,ii)}n>i a- stationary sequence with multivariate extremal index d'^{Ti, ...,Td), 
(ti, Trf) S Mq", and common distribution function Fx with unit Frechet marginals {x) ~ exp(— l/x), 

X > 0, and in the domain of attraction of G, i.e., Cp (11^", uV") — > Cg(wi, u,^). Consider 

an i.i.d. sequence of random vectors {Z„ = (Z„.i, Z„.(i)}n>i, also with unit Frechet marginals 
Fzj{x) = exp(— 1/x), a; > 0, and satisfying Fz 6 'D{H). Let a = [ai, ...,ad) be a vector of real 
positive constants. 

The pMAX model is defined by 

{Y„ = (K„.i, r„,d)}„>i = {(X„,i V \ ...,X„,d V Zll2')}n>u (9) 

with notation a V 6 = max(a, b). The letter "p" stands for the power transformation concerning {Z„}„>i. 
The EM4 model (Heffernan et ai, [5] 2007) is obtained by considering {X„}„>i an M4 process, aj = a, 
j ^ 1, d and Zn^i, Zn_d independent. 

At this point, we state some notation used along the paper. We will denote U the set of indexes j 
in D = {!,..., d} for which aj > 1. For any A C D, the vector (xi, ...,Xd)A denotes the sub-vector 
of (xi, ...,Xd) with indexes in A and, for any d.f. F, Fa denotes the marginal d.f. of F for sub- 
vectors with indexes in A. Also M„ and M„ will denote the componentwise maxima of, respectively, 
the pMAX sequence {Y„}„>i and the corresponding i.i.d. sequence {Y„}„>i. We will also use notation 
a Ab ~ min(a, b). 



2.1 The extremal index 

In this section we compute the marginal and the multivariate extremal indexes of a pMAX process. 
We will find the interesting and novel feature (when comparing with the EM4) that pMAX may have 
multivariate extremal index not constant equal to one while some marginal extremal indexes may be unit. 
First we need to compute the normalized levels form of pMAX. 

Proposition 2.1 A sequence of normalized levels {ui^^ = (ul^'j , 1 < j < d)}n>i for a pMAX process is 
such that, for each j E D, 

_,l/c 



h < 1 



.a, > 1 . 
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Dem. First observe that 



If aj < I, then 



l/a ■ j^l'^i 

1™ ^^^^ . ((-) ') " lim exp ( - I - = exp(-r,), 



'J 

whereas for a, > 1, we have 



n— J-oo \ \ 7"t / / n— >oo \ 7^ 



hm Fv ( — 1 = hm exp ( — r, ^ — t-\ ~ cxp(— r,). □ 

Now we compute the marginal extremal indexes of a pMAX process. 

Proposition 2.2 For each j G D, if aj > 1 then {Y„_j}n>i has extremal index coincident with the one 
of {X„.j}n>i, otherwise it will be unit. 

Dem. For the case aj < 1, we have normalized levels u^^'j = ("-/'''j)^^"^ fProposition l2.ip and hence 
Observe that 

lim^^oo i^x,^ ((;^) ') = hm„^oo exp - 4/^) = 1 

and thus 

lim„^oo P(Vr=i < iri/T,)^-^) = e-^^xo = 1. 

Therefore, 

lim„_^oo P{Mn,j < Un'j) = l™n^oo ■P(Vr=l ^ "/'^j) = S""^' ' 

If ttj > 1, the normalized levels u'J^^^ = n/r, (Proposition 12. ip and hence 

lim„^„o P{M^^, < uj^') = hm„^oo P{\/7=i < (^/^jO'/^O-PCVLi < "/^j) = e"'^"^ - ^ 

In what follows we adopt the convention that ifC/ = 0orD — J7 = the quantities involving the 
respective sub-distributions are considered null. 

Proposition 2.3 The multivariate extremal index o/{Y„}„>i is given by 

_ 9]^{Ti,...,Td)u\ogCGu{e-^\-,e-^')u + \ogCHu_u{e-^\-,e-^')D-u 
^^''■■■'""'^ logCG„(e--S...,e--)a + logCff,_,(e--S...,e--)z5_a ' ^ ' 



where 



9]f{Ti,...,Ta)u ^ , lim ^^(n, r^). 

TjlO, ]£D — U 



Dem. We have, successively, 

p({vr=i < (^)"^ J- e c^}, {vr=i < e ^ - 1^}) 
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Since the second terms in each of the two factors converge to zero, we have 

lim P(M„ < ul")) = hm F(V:Li X,,, < e U)P{\Jti < e D - U) 

f{ri,...,Ta)u , 

On the other hand, 



Um P(M„ < ul")) = hm P"(r„.j < j G D) 

= hm P"(^nj < ^, i e t/)P"(^n.J <^,J<^D-U) (11) 

Now the result fohows from relation ([7]) . □ 

Observe that, from expression (|TU)) . we obtain the result of Proposition 12.21 for the univariate 9j , 
JG{1,...,4. 

We can also conclude that, if aj > 1 for all j G D then the multivariate extremal index of {Y„}„>i 
coincides with the multivariate extremal index of {X„}„>i. If aj < 1 for all j E D then the multivariate 
extremal index of {Y„}„>i is unit. In the other cases, the multivariate extremal index of {Y„}„>i de- 
pends on the multivariate extremal index of the sequence {(Y(7)„}„>i of the sub- vectors with components 
having indexes in U. 

We shall see some particular cases: 

(i) If CG{ui,...,Ud) = Ch{ui, ...,Ud) = Aj=i"ii then 

(Tl,...,Td) = -— . 

V jeu 'j ^ V jeD-u '3 

(u) If Cg(ui, = Aj=i"j and Cniiti, ■■■,itd) = 11^=1%) then 

6l^(ri, Td)u V^g^ Tj- + T,jeD~u 



(n, ■■■,Td) 



(in) If Cg(wi, .■.,Ud) = Ch{ui, ...,Ud) = {Y,f=i{~ logUj)^/'^)!^), for some < ^ < 1, then 

, _ ^^(n, Tdhi r^^')' + ( E.eo-u 



(I]je(7'^j ) + (SjGD-C/'O' ) 

(iv) If aj = a, j € D, then 6^ {ti, r^) = 1 whenever a < 1 and 9^ {ti, r^) = ^■'^(ti, Td) if a > 1. 
This corresponds to the case of EM4 model (see Heffernan et ai, [5], 2007). 

2.2 Extremal coefRcients of the limiting MEV 

When computing the multivariate extremal index of the pMAX, we have already seen that if Fx G Z^(G) 
and Fz G V{H) then 



Therefore, for {xi, ...,Xd) E K'^, 

P({n-iM„,j < xj,j e [/}, {n-i/".M„,, <x,,jeD- U)) 

We shall denote 

nxi,...,Xd) = Gc/''"(^/"^-^/"^'^(xi,...,Xd)a-ffz3-c/(xr,...,x^^)z3_c/. (12) 

We have also seen in ([TT|) that i^Y G 2^(1^) where we have 

V{xi,...,xd) = Guixi, ...,Xd)uHD-uix"\ ...,x'^'')d-u- (13) 

The MEV distributions V and F do not have identicahy distributed marginals. Indeed, based on the 
proof of Proposition 12.21 we may see that, for aj < 1, 

P{n-^/°'^Mn.n < X.) — > e~^^ °' 

and for > 1, 

P{n-^M„.j < Xj) — > e-"'''^'. 

n— ^cxD 

Observe that we can now obtain these univariate extreme limiting distributions from the MEV ones, V 
and V. 

Since MEV distributions G and H have each identically distributed marginals, the definition of the 
extremal coefficient e of Tiago de Oliveira ([10], 1962/63) and Smith ([H], 1990) is pacific within this 
case. More precisely, for all x > 0, we have 

G{x,...,x) = {e-'"'')"' and H{x,...,x) = {e""'')'" 

and thus 

EG = - logG(l, 1) and en = - logif(l, 1). 

When the marginals of a MEV are not equally distributed, the previous definition needs some modifica- 
tion. A natural way is to consider an analogous definition based on MEV copulas. More precisely, the 
extremal coefficient of a MEV distribution V is the constant cy such that 

Cv{u,...,u)^u'^,ue [0,1]- (14) 
If V has identically distributed marginals, both definitions are equivalent. 

We shall calculate the extremal coefficients of V and V according to definition in (|T?)) . 

Proposition 2.4 The extremal coefficients of the MEV distributions, V in U^) andV in U3\}. are given 
by, respectively. 



and 

ev = 



Dem. Considering the variables change, Uj = V{xj) or Uj = Vj{xj), we obtain 

Cy{ui, ...,Ud) = Cgu{ui, ...,Ud)u CHo-ui^l^ ■■■,Ud)D-U, 

which leads immediately to the first assertion, and 

By the multivariate extremal index homogeneity property in ([8]), we have 

Cy (7., . . . , = C^;^ ^ ' • • • ' ^ ' " (yl/^S . . . , ^.1/'^- )a Cff,_ „ (u, . . . , , 

If we use MEV's spectral representation of G (Haan and Resnick, [6] 1977), then 



jeu 



u 



\l ^ dWG(wu...,W2) 



where Wq is a finite measure over the unit sphere 5*^ of R'J_. Now just observe that 
^v=0^{j^,-,j-Ju [ y J- dWG{wi,...,W2) + eHn-u- n 

Let us apply the result in the case where {X„}„>i is an M4 process. We have 

oo oo d 



CG{u^,...,ud)=x{ n A 

l — \ k— — oo j — 1 



where aikj. I > 1, —oo < k < oo, j E D, are the model constants and 

9 (Tl,...,Trf) = — ^ 



1=1 Z^fc=-oo y j=iU] aikj) 
Therefore, any choice of H will lead us to 

oo oo oo oo 

i=ik=-oojeu 1=1 k=-oo jeu 

where 

cxj oo 



A;— — oo 

If in particular we choose H as in Heffernan et al. 2007), we have Cy = ey = eh = d \i a < 1 and, 
for a > 1, 

oo oo oo oo 

^v^Y.H y and = ^ V V 

i=i fc=-oo jeu i=ik=-ooieu ^ 
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2.3 Tail (in) dependence coefficients 

In this section we characterize the tail (in)dependence of a pMAX model {Y„}„>i in We will com- 
pute the lag-r TDC, A^y (Y) and the lag-r Ledford and Tawn coefEcient, rj^yiY)^ by assuming that the 
respective exist for the underlying process {X„}„>i. 

For calculations, we shall need the pMAX bivariate d.f. 

PiYl.j < X],Yi+r.j' < X]>) = PiXi^j < Xj,Zij < x"' ,Xi+r.j' < Xj',Zi+r,j' < ) 



Proposition 2.5 // the pMAX process {Y„}„>i in (0) is such that A^y (X) exists for the underlying 
process {X„}„>i, then 

,aj<l 



Ag(X) 



Dem. We have 

P{Yi+rj' > x\Yij > x) 



1 - PjYl.j < X)- P{Yl+r.r < + P(yi,j- < X, Yi+r,j' < X) 



^ P{Xi+r,j' < x)e-"'°^' + {P{Xi+r.j' <X)~ P(Xi,, > X, Xi+r.y < x)) C 

1 - P{Xij < a;)e-^"°' 

-1 ~c., 1 _ e-^"°" -c- -c, 1 - p"^" 

= +P{Xi+r,r <x\Xi^,>x)e-- ' ' 



Now just observe that 

1 , aj <1 
1/2 

, aj > 1 
and 



lim — — = 0/2 , a, ^ 1 

X-¥00 \ _ p-X-'^-X J ' 



1 -1 I , a,- < 1 



lim : r^— ^ = <! 1/2 , a, = 1 

1 , Qfj > 1. □ 

In the particular case of the EM4 in Heffernan et al. ([5], 2007), we have aj = a, j = 1, ...,d, leading 

to 

, a < 1 

aS(Y) = <' U^-S) ,a = l 



2-S , a > 1. 



where 6 = X]fe°=-oo '^zfcj ^ '^ikj' (^^^ Result 9 in Heffernan et ai, jS] 2007; observe that a small 

correction is needed for the case a = 1). 
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Proposition 2.6 // the pMAX process {Y„}„>i in ^ is such that 77^^,'' (X) exists for the underlying 
process {X„}„>i, then 



max 



Dem. By hypothesis, 



P{Xij > x,Xi+r,j' > - L(a;)x"^/"^^'^^^ (15) 



for some slowly varying function L. 
We have 

P(ri J > X, Yi+rj' >x) = l~ e"^"'-^""' - e-^"'-^"°'' + P(Xi ^ < x, Xi+r,]' < x)e^^"°' 
where 

P(Xij < X,Xi+r,j' <X)= 26""=^' - 1 + P(Xi^j > X,Xi+r,j' > x). 

Observe now that, for any positive real constants a, b and large enough x, 

-x-'^-x-" 1 1 1 1 1 1 

e ~ 1 \ \ 1 

.j,a j.b 2x2" 2x2'' 2x°+''' 

Therefore, after some calculations, we obtain 



which implies the assertion. □ 

In the EM4 case, the process {X„}„>i is an M4 and thus ri^y{X.) = 1 and aj ~ a, for all j ~ 1, d. 
Hence, we have 



%'m = < V2' 

t 1 , a > 1 . 

already derived in Heffernan et al. ([5] 2007; Result 10). 

Observe that, in a pMAX process, aj < 1 leads to tail independent random pairs (Yj-,Yj') for all 
j' € D, whereas for aj >1, random pairs {Yj,Yj') will be tail dependent (independent) wether {Xj,Xjr) 
are tail dependent (independent). This is one important advantageous of the pMAX: it is suitable for 
multivariate time series data for which some pairs of components may present tail dependence while 
others may not. 

Example 1 Consider {X„}„>i = {{Xn^i, Xn^2, Xn,3)}n>i 0, i- dimensional sequence with X„.i = C/„, 
^n,2 = Wn and X„.3 = Wn+i, whcrc sequences U = {C/„}„>i and W = {tV„}n>i o,re such that: 

(i) U is i.i.d. with common d.f. Fij ; 

(ii) Fjj g 2?(Gi,{a„ > 0},{6„}), i.e., Fjj belongs to domain of attraction of a non degenerate extreme 
values distribution Gi with normalizing constants {a„ > 0} and {6„}; 
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(Hi) U is independent of a chain J = {Jn}n>i of independent Bernoulli(l/2) r.v.'s; 
(Hi) Wn = ?/«l{,/„=o} + ?7„+il{,/„=i}, n > 1. 

Observe that U and W have the same common d.f. and thus normalized levels u'n'^ for U are also for 
W. By hypothesis, F'^^ . {onX + 6„) — >■ Gi{x), j = 1, 2, 3, and we have 

1 ^ 

i^x„(xi,X2,X3) = -(2F[/(xi Aa;2)Fc/(a;3) +Fc/(xi)Fc/(a;2 Axa) + n^t'(2^j) 

i=i 

Now observe that 

n{l - ^x„ {anXi + bn, a„X2 + 6„, a„a;3 + 6„)) 

= f (4 - 2i^(7(a„(a;i A X2) + 6„)F(7(a„X3 + 5„) - Fu{anXi + b„)Fuian{x2 A 2:3) + 6„) 

= f (2(1 - Fu{an{xi A a;2) + 5„) + Fu{an{xi A 0:2) + &„)(1 - i^c/(a„X3 + 5„))) 

1 - Fu {unXi + bn) + Fu{a )(l-Fu{an{x2AX3)+bn)) (16) 

1 - F[/(a„a-i + bn) + Fu{anXi + 6„)(1 - Fu{anX2 + bn)) 

+Fu{anXi + bn)Fu{anX2 + 6«)(1 - i^t7(a„X3 + fo„))^ 

i( - logG'i(xi A 2:2) - logGi(a:;3)) 

+ i( -logGi(a-i) -logGi(a;2 AX3)) 

i ( - log Gi (xi ) - log Gi (a;2) - log Gi (0:3)) . 

T/iere/ore, 

-Px„ (ana;i + 6„, a„2:2 + 6„, a„a;3 + 6„) 
^ G(xi, X2, X3) = G}/'(xi A X2)G\'\x2 A .T3)G}/'(xi)G;/"(x2)G?/'(x3), 

n— >oo 

/rom which we derive the copula 

Cg{ui,U2, U3) = {ui A U2)^^^(w2 a U3)^/'^uY^ul^^ul^^. 

Next we compute the multivariate extremal index o/{X„}„>i, which is a 2-dependent sequence and thus 
satisfies a multivariate version o/D'''^'(u'-'^') condition of Chernick et al. (J^, 1991), 

lim ^(Xi ^ ''n\^2 < u*:', X3 < ui:\x, ^ uir^) = o, 

n— >-oo ^ — ^ 

/or any sequence {fc„}„>i suc/i f/iaf, fc„ — > 00, knan,i„ 0, fc„/„/n — > 0, where {an.z}n>i is sequence 
of strong-mixing coefficients and {ln}n>i is such that a„_;„ ^ 0. Therefore, the multivariate extremal 
index is given by the following limit, whenever it exists: 



71— ^00 



hm 7;^- =6* (Ti,r2,T3). 



Observe that the denominator can be derived from Hid]) , leading to 

lim nP(Xi ^ u^;)) = i(2(ri V ra) + (t2 V rg) + 2ti + r2 + 8x3). 



11 



Table 1: The lag-r TDC coefficients of sequence {X„}„>i. 





j'=l 


j'=2 


.i'=3 


j=l 


if r > 1 


1/2 if r = 
if r > 1 


if r > 


j=2 


1/2 if r = 1 


1/4 if r = 1 


1/4 if r = 




if r > 2 


if r > 2 


if r > 1 


j=3 


1/2 if r = 1,2 


1/2 if r = 1 


1/4 if r = 1 




if r > 3 


1/4 if r = 2 
if r > 3 


if r > 2 



For the numerator, after some lengthy but simple calculations we obtain 
lim nP(Xi ^ uM,X2 < uM,X3 < ui^^) = 

n— >-OD 

and thus we have 



4ri 




, Tl 


> 


T2 


> 


T3 


Tl 


> 


T3 


> 


T2 


Tl 4 


' 3t3 




> 


n 


> 


T2 




> 


T2 


> 


Tl 


Tl 4 


- 3t2 


, T2 


> 




> 


T3 


T2 


> 


T3 


> 


Tl 



e''(Tl,T2,T3) = 



as well 



^^2(Tl,T2) 



4ri 



4Ti+2r2 
2ti+4t2 



4ti 



4ti+2t2+3t3 
•ri+3T3 

4Ti+r2+4T3 

4ri 

4ti+T2+4t3 
T1+3t2 

2Ti+4r2+3r3 

Ti+3t3 

2Ti+3r2+4r3 



Tl > T2 
T2 > Tl 



Tl > T2 > T3 

T3 > Tl > T2 

Tl > T3 > T2 

T2 > Tl > T3,T2 > T3 > Tl 

T3 > T2 > Tl. 



^M(n,T3) 



T1+3t3 

4ri+4T3 

4ti 
4ri+4T3 



, T3 > Tl 
, Tl > T3 



3r3 



^2''3(^2,T3) 



3t2+4t3 

3r2 
4r2+3T3 



, T3 > T2 
, T2 > T3, 



and univariate extremal indexes, 0^ = 1 and 0^ = Of- = 3/4. 

Applying the lag-r TDC definition in simple calculations lead us to the results summarized in 
Table[l\ Obviously rj^-^-jCX.) ~ 1 whenever Ajy(X) is positive and, since we obtain P(Xi+,.j' > x,Xij > 



x) = P{Xij > x)^ in all cases of a null TDC, we have ?/,y^'(X) = 1/2 in those cases 



Now if we consider a pMAX process 



with ai = 3/2, a2 = I and = 2/3, then by applying the Provosition \2.3[ we have 9^ = 9^ = 3/4, 



.1/02 



9Y 

-^2,3 



(t2,T3) 



|t2+T3 
T2+r3 



and 9^{ti,T2,t^) 



2 (ti ,'''2 ) (ti Vr2 ) +ri +T2 +2t3 

(ti VT2)+Ti+T2+2r3 



4r^ +3ti T2 +4ti T3+2T2 T3 
4r +t| +4ti T2 +4ti T3 +2t2 T3 

2r j- +7t| + 7ri T2 +4ri r3 +8t2 T3 
2t'1 +8T|+8ri T2 +4ri r3 +8T2 T3 



T/ie pM^Z TOC's are A^'^.?(Y) = aJ'^-!(X), A^^-!(Y) = 5A^^-^(X) and A^^-!(Y) = 0, /or all r e No and 
j' = 1,2,3. In what concerns the lag-r Ledford and Tawn coefficients, we have if^-^ (Y) — 77^2' (Y) = 1/2 



Tl > T2 
T2 > Tl . 
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if j = 1,2 and r > j, -nja (Y) = 3/5 if j = 1,2 and r > j - I, r]^y(Y) = 2/5 if j = 1,2 and r > 3 and 
r,il\Y) = l/2 tfr>2. 

3 Estimation: some notes 

In this section we shall give some guide marks on the estimation of the marginal parameters aj, j G D, 
and the tail dependence coefficients within a pMAX process. For the estimation of the multivariate ex- 
tremal index, we refer the manuscript of Smith and Weissman ([12], 2006) and the paper of Robert (|15j. 
2008). 

We start by remarking that aj, j G D, are the key parameters of the pMAX process dictating the 
type of dependence and thus allowing to sometimes deduce the value of some coefficients. For instance, 
an estimate of aj less than one means that the random pairs {Yn_j, Yn+r,j'), for all n,r £ N and j' G D, 
are tail independent leading to a respective null lag-?- TDC, Ajy(Y). In addition, estimates of aj less 
than one for all j & D will imply an unit multivariate extremal index. So, as a first step, we suggest to 
marginally apply a test of tail independence (see e.g, Zhang [5T] 2008 and references therein). 

In order to estimate the value of aj, observe that it corresponds to the tail index of the marginal 
process {Ynj}n>i whenever aj < 1. There are several known tail index estimators in literature such as. 
Hill, Pickands, maximum likelihood estimator, moments, generalized weighted moments, among others. 
Their properties have been derived in an i.i.d. framework, but there are some studies considering a 
stationary context (see, for instance, Rootzen et al. [Tn| 1990, Resnick and Starica |131 [T^ 1995/1998, 
Drees [21 2003). The Hill estimator is the most used in the context of a Frechet domain of attraction. 
In particular, under a strong-mixing dependence structure, it is consistent and asymptotically normal 
(Rootzen et al. jT6j ) . 

For a given sequence {^n}n>i, a strong-mixing condition means that ^„ and ^„+s are approximately 
independent for large s. More formally, it will hold whenever 

\P{AnB) ~ P{A)P{B)\ <a{s) 

s—^oo 

for A G ...,^k), B G T{S,k+s,£.k+s+i,--) and fc,s > 1, where ...,£,k, —) denotes the a-algebra 

generated by {^i, ...}. This is a mild condition to assume as many stationary processes, including M4 
and EM4, do satisfy strong-mixing. In our case it will be easily stated for the marginal process {l'„j }„>i 
by just assuming that it holds for the underlying process {X„.j}„>i (Bradley ([T], 2005; Theorem 5.2). 

We can also find in literature estimation procedures of the TDC and of the Ledford and Tawn coeffi- 
cient. However, the properties of consistency and asymptotic normality are mainly derived under an i.i.d. 
assumption, i.e., considering i.i.d. random pairs, which is not an interesting case when we intend to esti- 
mate the respective lag-r versions, A^y (Y) and tI^^,{Y). A brief note on the estimation of the TDC, when 
we drop the independence assumption, can be seen in Ferreira and Ferreira ([S], 2012). In what concerns 
coefficient r{-^^,{Y) defined in (0]), observe that it can be estimated as the tail index of min(Fi.j, Yi+rj')- 
By assuming that the underlying sequence {X„}„>i satisfies the strong-mixing property we also derive 
a strong-mixing structure for sequence {m\n{Yn,j,Yn+r.j')}n>i (see Lemmas lA.ll and IA.2I in Appendix). 
Now the considerations above concerning the Hill estimator will also hold in this case. 

A Appendix 

Lema A.l Let {Y„}„>i he a d-dimensional sequence satisfying the strong-mixing condition with a se- 
quence of dependence coefficients {a(s)}s>i. Then, 

(i) {(yrij, ^n-Hrj"')}n>i Satisfies the strong-mixing condition with any sequence of dependence coeffi- 
cients {6(s)}s>i such that b{s) = a{s — r), s > r. 

(ii) {(Ynj A y„+rj"')}n>i Satisfies the strong-mixing condition with {6(s)}s>i o-s in (i). 
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Dem. 

(i) By hypothesis, V.4 e T{Yi, Yp), B e J'(Yp+s, Yp+,,+i, ...), s > 1, we have 

\P{AnB)- P{A)P{B)\ <a{s) 0. 

s— >-oo 

Consider r > 1 fixed, C £ ^'((Yij, Yi+^j/), (ypj, Fp+^j/)), £> e T{{Yp+s.j ^Yp+r+sj'), ■■■) and 
s > r. We have C E J-{Yi, Yp+r), -D G J'(Yp+s, Yp+^+i, ...) and hence 

\P{CnD)-P{C)P{D)\<a{s^r) ^ 0. □ 

(ii) Just observe that 

for any family of indexes /. □ 

Lema A. 2 Let {X„}„>i be a d- dimensional sequence satisfying the strong-mixing condition and {Z„}„>i 
a d-dimensional i.i.d. sequence independent of the previous. Then {Y„ = X„VZ„}„>i satisfies the strong- 
mixing condition. 

Dem. Consider sequence {U„ = {X„ i,...,Xn,dTZn^i,...,Zn ,d)}n>i and let A e J'(Ui, Up), B e 
J^(Up_|-s, Up_)_s+i, ...). Since we can write 

\P{A riB)- P{A)P{B)\ = \P{A' n B')P{A" n B") - P{A')P{A")P{B')P{B")\, 

with A' e ^(Xi,...,Xp), B' e .F(Xp+,„Xp+,+i,...), A" e J-(Zi,...,Zp), B" e ^(Zp+„ 
have that 

|P(A n B) - P{A)P{B)\ < \P{A' n B') - P{A')P{B')\P{A" n B") 

+ n B") - P(A")P(B")I^'(^')^(5') 
< ax„(s) ^ 0. 

Therefore sequence {U„}„>i satisfies the strong-mixing property. Now the result follows from the same 
argument used in Lemma lA.ll fii) for the maximum operator (V). 
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